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This Letter reports on how the interfaces in the (l+l)-dimensional Kardar-Parisi-Zhang (KPZ) 
class undergo, in course of time, a transition from the flat, growing regime to the stationary one. 
Simulations of the polynuclear growth model and experiments on turbulent liquid crystal reveal 
universal functions of the KPZ class governing this transition, which connect the distribution and 
correlation functions for the growing and stationary regimes. This in particular shows how interfaces 
realized in experiments and simulations actually approach the stationary regime, which is never 
attained unless a stationary interface is artificially given as an initial condition. 
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Aside from their ubiquity in nature, surface growth 
phenomena constitute an important situation of statis- 
tical mechanics out of equilibrium, where scale invari- 
ance and universal scaling laws arise generically [T]. At 
the heart of such growth processes is the Kardar-Parisi- 
Zhang (KPZ) equation [5] and the corresponding univer- 
sality class [H describing the simplest case without 
any conservation laws and long-range interactions. For 
one-dimensional interfaces, the KPZ equation reads 

3,, , d 2 h A / dh\ 2 /— , 

where h(x, t) denotes the fluctuating height profile and 
r](x,t) white Gaussian noise with (rj(x,t)) = and 
(rj(x,t)r}(x',tf)) = 5(x - x')5(t - t'). The values of the 
scaling exponents are exactly known in this case [THS]: 
the height fluctuation dh = h — (h) grows with time as 
Sh ~ t 1 / 3, and the correlation length £ as £ ~ t 2 ! 3 . Specif- 
ically, h is described by an appropriately rescaled random 
variable x(a/,i) as 

h(x 1 t)^v oc t+(Tt) 1 ^ X {x',t) (2) 

with a rescaled coordinate x' = (Ax/2)(Tt)~ 2 / 3 and con- 
stant parameters A = v/2D, T = A 2 \\\/2, and v^. The 
KPZ-class exponents have indeed been reported in var- 
ious models and theoretical situations [IHl]j as well as, 
recently, by a growing number of experiments [5l410|. 

Studies on the (l+l)-dimcnsional KPZ class entered 
an unprecedented stage in 2000, when Johansson [IT] and 
others [4 rigorously derived asymptotic distributions of 
the height fluctuations for a few models. Among oth- 
ers, it has brought about two outstanding outcomes: (i) 
The KPZ class splits into a few subclasses according to 
the global geometry of the interfaces, or, equivalcntly, 
to the initial condition. These subclasses are character- 
ized by different distribution and correlation functions, 
whereas they share the same scaling exponents, (ii) An 
unexpected link to random matrix theory has been re- 
vealed. In particular, the asymptotic distribution of \ 
for the flat and curved interfaces is given by the largest- 
eigenvalue distribution, called the Tracy- Widom (TW) 



distribution [T2], for the Gaussian orthogonal and uni- 
tary ensembles (GOE and GUE), respectively |T3 . The 
stationary interfaces also form a distinct subclass. To 
study it analytically, one usually sets the initial condi- 
tion h(x, 0) to be a stationary interface, which is sim- 
ply the one-dimensional Brownian motion for the KPZ 
equation [JJ. The height difference h(x,t) — h(x,0) then 
grows as Eq. (J2| and \ obeys the so-called Fq distribu- 
tion, as proved for the polynuclear growth (PNG) model 
[131 114] , for the totally asymmetric simple exclusion pro- 
cess P21 US], and, very recently, for the KPZ equation 
[17] , The two-point correlation function being exactly 
derived as well [T5HT9] , this subclass is now firmly estab- 
lished like the ones for the flat and curved interfaces. 

Such a stationary regime is, however, never attained 
within a finite time in an infinitely large system, unless a 
stationary interface is artificially given as an initial condi- 
tion. This is readily seen by recalling that the correlation 
length grows as £ ~ t 2 / 3 , whereas it is infinite for the 
stationary interfaces. Therefore, in practical situations 
starting from a smooth or uncontrolled initial profile, one 
needs to elucidate how the interfaces approach the sta- 
tionary regime in course of time. This is achieved by the 
present Letter. Simulations of the PNG model show that 
the height difference of the flat interfaces exhibits a tran- 
sition from the flat, growing regime to the stationary one. 
We find scaling functions describing this crossover and 
determine their functional forms, which smoothly con- 
nect the GOE-TW and F distributions. We also study 
the two-point correlation function and show how those 
for the two subclasses interplay at finite times. These re- 
sults are quantitatively reproduced by experimental data 
on turbulent liquid crystal pj [5] , indicating that they are 
universal characteristics of the KPZ-class interfaces. 

First we numerically study the PNG model. Starting 
from a flat substrate h(x, 0) = 0, an interface experiences 
random nucleation events at a uniform rate. On each nu- 
cleation, the local height h{x, t) increases by one, produc- 
ing a plateau that expands laterally at constant speed. 
When two plateaux encounter, they simply coalesce. The 
interface grows thereby upward on average. We numer- 
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FIG. 1. (color online). Crossover in the one-point distribution for the PNG model. (a,b) First- to fourth-order cumulants 
(Ag n ) c against At (a) and At/to (b), for to = 0,0.1,1,10,100,1000,4000,7000 (increasing as the arrows indicate). The top 
(dotted) and bottom (dashed) horizontal lines indicate the corresponding values for the Fq and GOE-TW distributions, (xo)c 
and (x")c, respectively. The black solid lines in (b) show fitting to the collapsed curves (see text). The insets in (b) show the 
skewness and the kurtosis. (c,d) Asymptotic behavior of the data in (b) for small and large At/to- The data for to = 0.1 are 
omitted because of the strong finite-time effect, (e) Distribution of Aq for given pairs of to and At/to. 



ically implement space-time representation used for the 
analytic derivation of the distribution function [T31 [T3] , 
with the nucleation rate 2 and the plateau expansion 
speed 1. This choice of the parameters corresponds to 
Voo = 2, A = 2, and T = 1. For the simulations, we 
impose the periodic boundary condition with system size 
10 3 and realize 10 4 independent interfaces up to time 10 4 . 

The quantity of interest is the height difference 
Ah(x, At, to) = h(x, tg + At) — h(x, to), rescaled here as 

Aq(x,At,t ) = A ( r^T/f • (3) 

By construction, Aq A xr for ~^ and then At — > oo, 
while Aq A xo for to — >• oo and then At — > oo, where xi 
and xo are random variables obeying the GOE-TW and 
Fo distributions, respectively, with the factor 2~ 2 / 3 mul- 
tiplied with the usual definition for the former [T31 [H] . 
Figure [ija) shows the first- to fourth-order cumulants 
of Aq, (Aq n ) c , as functions of At for different to, dis- 
played with the corresponding values for the GOE-TW 
and Fo distributions. The cumulants agree with those for 
the GOE-TW distribution as At tends to infinity, while 
they indicate the values of the Fo distribution for large 
to and small enough At. The transition curves are found 



to collapse very well when At is scaled by to [Fig. [TJb)] , 
except for the early stage with too small to and At. In 
particular, for to — > oo, the cumulants converge to a sin- 
gle set of functions, (Aq n ) c — > AQ n (At/to), satisfying 
AQ„ {t) -> (xi) c for t -> oo and AQ„(t) (xo) c for 
r — > 0. One can indeed draw the functions AQ n (r) by 
making histograms for (Aq") c at each At/to with vary- 
ing to and fitting their modes by, e.g., spline functions, as 
shown by the black solid lines in Fig. [ljb). Theoretical 
expressions of AQ n (r) are not obtained, because they in- 
volve the time correlation which has not been solved yet. 
Asymptotically, the data suggest (xo)c — AQ\{t) ~ t 2 / 3 , 
(xjj)c - AQ 2 (t) ~ r 1 / 2 for small r and AQi(r) - (vi) c ~ 
T-V3, AQ 2 (t) - (xl)c - r- 2 / 3 for large r [Fig. |ljc,d)], 
while one needs better statistical accuracy to determine 
asymptotic laws for n > 3. In between the two limits, 
the transitions occur earlier for larger n (< 4), leading to 
interesting undershoot in the skewness {Aq 3 ) c llAq 2 )^ 2 
and the kurtosis (A g 4 ) c /(Aq 2 ) 2 [insets of Fig. [gb)]. Fi- 
nally, this crossover and scaling in terms of At /to can also 
be checked directly in the distribution; Fig. [TJe) shows 
that the probability density functions of Aq overlap for 
fixed At I to and that they shift from the Fo to the GOE- 
TW distributions as At /to is increased. 
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Now we turn our attention to the two-point correlation 
function, denned here by 

C{1, At, t ) = ( [6h(x + l,t + At) - Sh(x, t )] 2 ) (4) 

with 5h(x,t) = h(x,t) — (h(x,t)). If one takes the 
stationary limit to —¥ oo and then considers large At, 
one has C'((,At,t ) = {TAt)- 2 / 3 C(l, At,t ) ~ g(() 
with C = (Al/2)(TAt)- 2 / 3 , where g(() is the exact so- 
lution for the rescaled stationary correlation [T71 1TB] . 
This is tested in Fig. [2]ja) with finite to and At, where 
AC"(C, At, t a ) = C*'(C, At, t ) - C*'(0, At, t ) is compared 
with g(£) — g(0) in the main panel. First we note that 
the data for fixed At /to and different to overlap with 
each other, confirming that At /to is the only time scale 
that controls the dynamics. Now, we focus on the data 
with the smallest At/to we have, namely At/to = 0.006, 
shown by solid symbols in Fig.[2ja) (top data set). When 
the rescaled length £ is small enough, the data agree with 
the stationary correlation function, — 5(0), as ex- 
pected above. This agreement is confirmed without sub- 
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FIG. 2. (color online). Crossover in the correlation func- 
tion for the PNG model, (a) AC'((,At,t ) against ( for 
At/to = 0.006,0.1,1,10,100 (different sets of data; increas- 
ing from top to bottom) and t = 1, 10, 100, 1000, 4000, 7000 
(different colors and symbols; increasing from right to left). 
The dashed and dotted lines indicate the stationary correla- 
tion g(()-g(Q) and the Airyi correlation 2(At/t )~ 2/3 [pi(0)- 
5i((Ai/t ) 2/3 C)], respectively, with At/t = 0.006 for the lat- 
ter [see Eq. (pjl]. The inset shows the data for At /to = 0.006 
without subtraction of C"(0,At,to), compared with g(Q. 
(b,c) Asymptotics of AC'((, At, t ) for £ ->• 00 (b) and ( ->• 
(c). Blue circles are the numerical estimates, obtained from 
the data with to = 10 for (b) and with various to for (c). 
The gray dots in (b) show the values of the right-hand side of 
Eq. B, where (Ag 2 ) c for t = 10 and 100 is used as (^(in- 



fraction of C"(0, At, to) as shown in the inset. In con- 
trast, for large C, AC"(C, At,t ) deviates from it and 
converges to a constant. This long- length behavior is 
actually governed by the spatial correlation of the flat 
interfaces, namely the Airyi correlation gi(-), defined by 
gi (v) = (Ai(u + v)Ai{u)) — (Ai(u)) 2 with the Airyi pro- 
cess Ai(u) [2T)H2"2"] . To see this, we take At ->• in Eq. Q 
and obtain, for large t , AC"(C,0,t o ) = C"(C,0,t o ) ~ 
2(At/t )- 2 / 3 [ ff i(0)- 5 i((At/t ) 2 / 3 C)]. This function 
with At /to — 0.006 is indicated by the dotted line in 
Fig. [2^ a) and accounts for the data with large In short, 
when At /t is small enough, the data indicate 
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(C«Cc), 
(C»Cc), 

(5) 

where the crossover length £ c is defined by the intersec- 
tion of the two functions in Eq. When At/t is fur- 
ther decreased, the Airyi branch in Eq. ^ moves away 
as (At/t )~ 2 / 3 along both axes, leaving, asymptotically, 
only the stationary correlation g(() as expected. Alter- 
natively, if C and I are rescaled by t 2 / 3 instead of At 2 / 3 , 
what remains asymptotically is the Airyi correlation. For 
tiny but finite At /to, the crossover between the two limits 
is simply described by the function C ~ 2£. 

We then study how the correlation function varies for 
large At /to- The data series in Fig. [2^ a) for each At /to 
lower their values with increasing At/to. In the long- 
length limit ( — > 00, since (5h(x + 1, to + At)5h(x, to)) —> 
0, we have AC"(C, At,t ) ->• 2(rAt)~ 2 / 3 C t (At, t ) with 
C t (At,t ) = (Sh(x, t + At)Sh(x, t )), i.e., the time corre- 
lation function. Though it remains analytically unsolved, 
theoretically, its short-time behavior (At /to <C 1) is given 
by 



c t (At,t )~ (r 2 t t r )3(xi) c 



(x 



2 ) 

0/c 



2(X?) 



1 - 



(6) 



with t r = t + At [SI [H For At/t > 1, numerical 
[24] and experimental [8] studies showed C t (At,t ) ~ 
(r 2 t t r ) 1/3 i ;l (At/to) with F(r) - r" 1 . They indicate 



{(At/to)- 2 / 3 , (At/to «1), 
\ (At/to)- 4 / 3 , (At/to »1), 



(7) 



lim AC7'(C,At,t ) 

7" OO 



and correctly account for the data, as shown in Fig. |2| 
Further, since the second-order cumulant of the rescaled 
height difference, Q2(At/to), involves the two-point time 
correlation Ct(At,to), we also obtain for arbitrary At/to 

lim AC"(C,At,t ) 



= (xi) 



At/t c 



[At/to 



-Q 2 (At/t ). (8) 
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This is also confirmed as shown by gray dots in Fig. [2^b) . 

In contrast to the long-length limit, one cannot a 
priori predict how the short-length limit £ — > of 
AC{(,At,t ) depends on At/t . The data in Fig. &&) 
suggest AC {C, At, t ) ~ C 2 for any At/t Q . Figure Qc) 
shows that the coefficient of this quadratic term varies as 



hmAC'(C,AMo)C 2 =2 ^ 



,g"(0)/2 fa 1.085 
c(At/t )-^ 



{At /t « 1), 
{At /t » 1), 



(9) 



where g"(0) is the second derivative of g(Q at ( = 
and c is a constant. The value for At/t — > is 
due to C{(, At, to) — > g{C) in this limit. To examine 

oo, let us note | 



the other limit At /to 

(A/2)- 2 (rAt) 2 / 3 



At)§{x,t )), which 



is sim- 
ply time correlation in the slope of the interface. It 
is suggestive that the short- and long-length limits 
of AC (£, At, t ) are governed by the slope-slope and 
height-height time correlations, respectively, decaying 
with the same power in the rescaled units [Eqs. Q and 
djj])]. The results may also remind us of the space-like and 
time- like paths argued in the literature [25] [26], along 
which correlation is essentially governed by the purely 
spatial and temporal one, respectively, though precise 
relation to this concept is yet to be clarified. 

Finally, we test universality of the crossover reported 
so far, analyzing experimental data of growing interfaces 
in turbulent liquid crystal. While detailed descriptions 
of the experiments are given in Refs. H], in this se- 
ries of work the author and a coworker studied planar 
evolution of borders between two distinct regimes of spa- 
tiotemporal chaos, called the dynamic scattering modes 1 
and 2, in the electroconvection of nematic liquid crystal. 
The interfaces grow under high applied voltage, clearly 
exhibiting, besides the exponents, the distribution and 
correlation functions for the flat and curved KPZ-class 
interfaces 011]. Here, we employ the data for 1128 flat 
interfaces used in Ref. [5] and perform the crossover anal- 
yses developed in the present study. 

Figure[3]shows the results. The nth-order cumulants of 
the rescaled height difference Aq [Eq. with various to, 
which sufficiently fall apart as functions of At [see, e.g., 
inset of Fig. Pjja)], collapse reasonably well when plotted 
against At /to [Fig.|3ja)], despite rather strong finite-time 
effect for n > 2. The collapsed data are found asymptot- 
ically on top of the fitting curves obtained for the PNG 
model, Q n {At/to) (black solid lines). This implies that 
Q n {r) are universal functions of the KPZ class describing 
the crossover in question, and so is the distribution func- 
tion of Aq parametrized by At /to- The undershoot in 
the skewness is also confirmed experimentally [Fig.[3][b)], 
while it was not clearly identified for the kurtosis because 
of larger statistical error (not shown). Moreover, extrap- 
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FIG. 3. (color online). Crossover in the liquid-crystal exper- 
iment. (a,b) Cumulants (a) and skewness (b) against At/to 
with t = 2, 6, 10, 18, 30, 54, 60 s (from right to left). The inset 
shows (Aq) c against At. The dotted and dashed lines indicate 
the values for the Fo and GOE-TW distributions, respectively. 
The solid curves show the fitting to the PNG data obtained 
in Fig. |TJb) . (c) Rescaled correlation function AC' At, to) 
against £ for given to and At /to as indicated in the legend 
(At /to increases from top to bottom). The red dashed and 
black dotted lines indicate the stationary and Airyi correla- 
tion functions as described in Eq. |BJ, respectively, the latter 
being set with At /to = 0.006. The three black solid lines trace 
PNG data in Fig. [5|a) with corresponding values of At/to. 



olation of the finite-time corrections in the cumulants 
allows us to roughly estimate the time needed for direct 
observation of the Fo distribution, longer than 10 3 s here, 
which is unreachable in the current setup [S] ■ 

The results on the correlation function are also repro- 
duced experimentally [Fig. [3jc)] . The functional form 
is parametrized solely by At /to (see two data sets for 
At /to = 10 _1 overlapping with each other) and agrees 
very well with the one obtained for the PNG model (black 
solid lines). In particular, the crossover between the sta- 
tionary and Airyi correlation functions [Eq. (J5|] is clearly 
confirmed for small enough At /to (top yellow data set). 

In summary, we have studied the flat-stationary 
crossover in the KPZ class, which takes place gradually in 
time. Though crossover in space has been studied occa- 
sionally 14. , crossover in time is identified here for the first 
time to the knowledge of the author. Analyzing numeri- 
cal and experimental data, we have found and determined 
universal functions describing the cumulants and the two- 
point correlation during this crossover. These functions 
show multifaceted relations to the analytically unsolved 
time correlation, and hence may provide an important 
clue toward its solution. Besides this fundamental im- 
portance, our results also indicate how interfaces realized 
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in experiments and simulations approach the stationary 
regime, which is never attained without full control on 
the initial condition. 
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